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Abstract 

Collective transport through a multichannel disordered conductor in contact 
with charge-density- wave electrodes is theoretically investigated. The statis- 
tical distribution function of the threshold potential for charge-density wave 
sliding is calculated by random matrix theory. In the diffusive regime weak 

localization has a strong effect on the sliding motion. 
PACS numbers: 73.20.Fz, 71.45.Lr 
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Wave (Anderson) localization is an important concept of the quantum theory of electron 
transport. It is well established by now that a random disorder potential may cause bound 
states with an exponential decay governed by the localization length. When this localization 
length is large compared to other length scales in the system the corrections to semiclassi- 
cal transport are small and the term "weak localization" is appropriate. Weak localization 
is well understood for normal Fermi liquids. However, localization in heterostructures like 
superconducting/normal metal (S/N) or Josephson (S/N/S) junctions is still actively inves- 
tigated, experimentally and theoretically. For example, the electron-hole coherence induced 
in a normal metal by a superconducting reservoir significantly modifies the weak localization 
correction to the electron transport through the hybrid structure Q. 

A charge-density wave (CDW) is a collective state similar to BCS superconductors (S) 
0. For instance, CDW reservoirs induce a proximity effect in normal metals (N) analogous 
to that in N/S junctions when the interface is parallel to the CD wave fronts (normal to 
the chains). In such a geometry the phase coupling through a defect-free CDW/N/CDW 
junction supports a sliding CDW mode ||. A process analogous to Andreev scattering at 
an N/CDW interface has been predicted ||. Advances in thin film growth and patterning 
of CDW compounds []5J enables the study of CDW transport in various geometries on a 
mesoscopic scale Quantum oscillations have been identified in the magnetoresistance 

of the CDW compound NbSe3 in the presence of columnar defects ||. Electrical contacts 
which inject the current parallel to the CDW chains allow detailed studies of the subgap 
conductance ||. 

Here we report characteristic effects of weak localization on CDW sliding through a disor- 
dered conductor. In contrast to the extensive literature on semiconductor- or superconductor 
hybrid structures this topic has to the best of our knowledge not yet been investigated. 
The statistical properties of diffuse conductors including the weak localization correction 
to the average conductance are well described by random matrix theory |T0 |. We calculate 



the statistical properties of the pinning energy in the presence and absence of time-reversal 
symmetry by averaging over random transfer matrices. We find that the pinning strength 



is much more sensitive to the localization-enhanced backward scattering than the normal 
conductance. Weak localization is destroyed when time reversal-symmetry is broken by a 
magnetic field ||11|| , leading to a substantial negative magnetoresistance. 

We consider a disordered metal with length L and width W sandwiched between two 
CDW reservoirs consisting of N half-infinite chains (Fig. [I]). The CDW state is characterized 
by a complex order parameter A exp(z%), where the amplitude 2 A is the gap in the electronic 
spectrum and its phase x denotes the position of the CDW relative to the crystal lattice. 
The amplitude of the order parameter A is taken to be equal on all chains and constant up 
to the interfaces. In the limit of large correlation lengths £ = Hvf/A ^> L, the Coulomb 
interaction between the chains and the phase coherence over the region allows us to define a 
single CDW phase x f° r the whole structure. For simplicity the material parameters of the 
CDW and the normal metal are taken to be matched at the interface. This corresponds to 
a normal metallic island fabricated by ion damage || or selective abrasion |TI| of the CDW 
material. Mismatch between a CDW and e.g. a true 2D or 3D metal or imperfect interfaces 
cause additional scattering which, in principle, can be incorporated into our formalism. 

Many semi-conducting CDW materials are characterized by a nearly one-dimensional 
electron band structure. On the other hand, there are also CDW materials with substantial 
transverse dispersion, like NbSe3. For the latter systems and in the diffusive (L ^> I) and 
quasi-one-dimensional (L 3> W) limits an electron flux from one chain is distributed evenly 
among the outgoing chains. These two "anisotropic" and "isotropic" limiting cases (see Fig. 
|]) are treated below separately. Interpolation between these two extremes interpolation 
seems reasonable. 

In a Wannier site representation with respect to individual chains, equilibrium properties 
of CDWs can be described by the retarded and advanced quasi-classical Green functions 
(g a/3 )fj A (x, e), where a,f3 = {R, L} and i,j = {1,...,N} [fT3"| . The subscripts {R, L} refer 
to the right and left moving electrons on the linearized branches of the electronic spectrum, 
and i, j are the chain indices. The 2N x 2N matrix Green functions Left (-) and Right (+) 
of the scattering region are related by the boundary condition 
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g+ = g 



(1) 



where M is the 2 N x 2N transfer matrix of the disordered region which satisfies the condition 
M^EaMEa = 1 in order to ensure current conservation, and M* = SiMSi in the presence 
of time-reversal symmetry. Here Si and S 3 are defined as 



(2) 



Additional scattering at the interfaces mentioned above can be incorporated by substituting 
M — > MrMM/,, where Mrl are the transfer matrices of the right and left interfaces, 
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respectively. The transfer matrix can be decomposed as |JTT 
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^ I (l + A) 1 ^ A 1/2 ^ 7 
A 1 ^ (i + A)V2 

where a' 1 ' with i = {1, 2, 3, 4} are N x N unitary matrices and A is a real diagonal matrix 
with elements < {Ai r .., A^v} < oo, frequently referred to as eigenparameters of M. In the 
presence of time-reversal symmetry ot^* = ct^ and ct^* = cr 4 ) . 

Backscattering represented by the non-diagonal elements of M pins the CDW condensate, 
i.e. no current flows in linear response. By forcing the system out of the equilibrium position 
Xo a chemical potential difference between left and right side arises. The threshold potential 
fix at which the condensate starts to slide corresponds to the maximum value of the chemical 
potential difference as a function of the condensate position x- 



lPr/2 = max [//+(%) - fi-{x)] 



(4) 



— fi-] (x) can be identified as the periodic pinning potential. In terms of the retarded 
and advanced Green functions the chemical potentials fi± of left and right CDW's read 
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The pinning potential is expressed in terms of the transfer matrix elements via the Green 



functions g± , which can be obtained from their equation of motion J13J and the boundary 
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condition Eq. (p]). In the energy integral of Eq. fl5|), the continuum contribution to the 
energy integration dominates the bound states contribution when ep ^> A. In this limit 
the bound states on the island may be disregarded and we arrive via Eq. at the simple 
expression: 

A 



Tr(V — 



(6) 



where r' = a^^-^a^ 3 ^ and r = — oS^ jt^oc^ are the (N x N) reflection matrices of the 



disordered region and 7 1 is the dimensionless electron-phonon coupling constant |L4 
The threshold potential (||) is proportional to a factor containing CDW material parameters 
and a "form factor" depending on transfer matrix elements of the disordered region. In 
contrast to the reflection probability Trr^r, fi? is not invariant under a transformation r — > 
urv, where u, v are arbitrary unitary matrices. This is a direct consequence of the broken 
translational invariance of CDW's in the presence of defects. Naturally, the threshold is 
positive definite. By extending it symmetrically to negative threshold values the distribution 
function can be constructed from the averages of all even moments. 

The theory of random matrices is first employed for the isotropic limit described above. 
The ensemble average of a quantity /(M) of a collection of random conductors with length 
L is defined as 

(/> = J f(M)p L (X)d^(M) (7) 
where d/x(M) is the invariant or Haar measure and pz,(A) is the isotropic probability density 



n[. In terms of the parameterization in Eq. (H), the measure is explicitly given by 



2/3 

d/i(M) = Jp(\)d N \ [J (8) 

where the Jacobian 

■MA) = ni A »- A /> ( 9 ) 

a<b 

and dii(ay>) are the invariant measures of the unitary group. The symmetry index (3 = 1 
describes the orthogonal ensemble in the presence of time-reversal symmetry, and the unitary 



ensemble j3 = 2 implies broken time-reversal symmetry. The isotropic distribution function 
P(X,L) = pl(\)J(\) satisfies an iV-dimensional Fokker-Planck equation with the mean-free 
path I as single parameter JTTJ] . From the calculated average in Eq. ([?[) it is possible to 



extract the distribution V{\xt) of the threshold indirectly via 

(/)= / diPrfMVQjr). (10) 
Jo 

In the following () and () A denote the average over unitary matrices and eigenparameters, 
respectively. 

We start by considering the ensemble averages of the 2n-th powers of Eq. (|6|). All even 
moments (/iy 1 )^ can be generated from the quantity 

by a Taylor expansion in rj. It is convenient to write this as a Gaussian integral over the 
generating variables tp and ip* 

e r/Tr(r'-rt)Tr(r'+-r) = 1 f j V | 2 /^Tr(r't -r)+ V *Tr(r'- r t) _ ^\ 

In this form we can average r and r' separately over the unitary matrices, since the num- 
ber of normal and Hermitian conjugated matrices must be equal in order to give non-zero 
contributions to the average [ 15] . Equation flTT| ) may therefore be rewritten as 

n2 V 09) 
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The first order term in the summation (n = 1) can be evaluated exactly 



(13) 



TrrTrrt\ (/3) = (R)W , (14) 

where we used R = Trr^r =Tr A(l + A) -1 . This result agrees with the microscopic calcula- 



tions of Dorokhov |Ib| for the statistical distribution of the reflection matrix of a thin wire 
with Gaussian white noise potentials and inter-chain hopping. From this relation we learn 
that the threshold potential is determined by the (diagonal) reflection probabilities (I^mI 2 )^ 

6 



at each chain. The (3 dependent prefactor reflects the enhancement of backscattering by a 



factor 2, which is suppressed when time-reversal symmetry is broken [O]. It is instructive 



to compare this result with the reflectance of the disordered region (TrrtrV j = (R)f 



'A 

which differs only by the absence of the prefactor. The localization correction to the normal 
resistance recognized by its /3-dependence is called weak because it vanishes to lowest order 
in N. In contrast, the pinning due to CD modulations induced by the reservoirs is strongly 
affected. In the limit N ^ 1 the (squared) pinning strength is doubled (f*r) = 2 (/Jj,)^ 
by the weak-localization enhanced backscattering at each chain. 

For the n-th power ^ (Tt r^Tr ^> we found analytical expressions only for large N. In 



the Gaussian limit [17 



i 7 r Er(r / -rt)TV(r / t-r)\ ^ 



(15) 



P Vs + L 

to leading order O(N ). Here, the average over the eigenparameters A has been performed 
to leading order in N using the classical solution } = (Ns/ (s + l)) n with s = L/l ]10|] . 

Comparing the coefficients with powers of rj on both sides, we finally arrive at a Wigner-like 
distribution function 

rP(fPr) = ¥ge-<"*'te), (16) 

for < n T < oo, and = 2 (A/iV 7 ) 2 (s/s + 1) . The mean value and the standard 

deviation a are easily calculated as (/ir)^ = A/N"yJiTs/f3(s + 1) and a = \/3 ~ 1.732. By 



breaking time-reversal symmetry a magnetic field destroys the weak-localization enhanced 
backscattering, thus reducing the average threshold potential by a factor of l/y/2 ~ 0.71 
even though the normal resistance remains unchanged! 

We now examine the anisotropic limit in which interchain scattering is absent. Time- 
reversal symmetry cannot be broken by a magnetic field so we can restrict the discussion 
to the j3 = 1 ensemble. The distribution of the unitary matrices is not random but sharply 



centered around the unit matrix. Therefore we take the unitary matrices to be diagonal. 
For a weak impurity A < 1, located on a single chain N = 1 we obtain in the ballistic limit 
(s « 1) : 

V a {vt) = 2 e -/4/(/4) iV = 1, (17) 

(t4) 



with (/i|.) = (2A/7) 2 s. The distribution of pinning energies of a defect localized to a single 



chain satisfies a "half" Gaussian distribution, with an average (/z^) = (2A/7) ys/ir. In the 
large- N expansion the probability for zero threshold vanishes according to 

VaM = |f) e ~ /4/(/4) N h (18) 

where = jj (J^^ (j^pA for arbitrary lengths s. 

A negative magnetoresistance has been observed in NbSe3 in the presence of columnar 
defects, which is caused by the reduced threshold potential in high magnetic fields ||. In 
semiconductor CDW materials as Rb .3Mn0 3 or TaS 3 the effect has to our knowledge not 
been reported. These observations are consistent with our findings since the anisotropy of 
the normal state conductance of NbSe 3 is low, in contrast to the fully gapped, strongly 
anisotropic CDW materials. 

Our results also support previous work related to the periodicity of the Aharonov-Bohm 



oscillations |18fl . As shown here, the average threshold potential of a multichannel disor- 
dered conductor is sensitive to weak-localization trajectories. In the weak pinning limit, the 
effective threshold potential of many uncorrelated columnar defects as modelled by metallic 
islands is given by its averaged value. Therefore, we can argue that the latter is periodic 
in h/2e, since this component reflects the contribution from time- reversed paths. Finally, 
we point out that in the presence of many uncorrelated metallic islands the distribution of 
thresholds smears out the threshold singularity in the current- volt age characteristics, and 
determines the functional behavior of the non-linear conductance. 

We conclude by summarizing our results. We present expressions for the threshold 
potential for collective CDW current through a multichannel disordered conductor in terms 

8 



of its reflection matrices. By means of random matrix theory we calculate the statistical 
distribution of the pinning strength with and without time-reversal symmetry. We find that 
localization which is weak on the normal conductance has a strong effect on the threshold 
potential. The breaking of time-reversal symmetry by a magnetic field reduces the threshold 
field, depending on the isotropy of the disordered conductor. The magnitude of the reduction 
ranges from zero in the anisotropic limit to l/y/2 in the isotropic limit. 
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We acknowledge support from the NEDO joint research program (NTDP-98). 



9 



REFERENCES 



[1] Transport Phenomena in Mesoscopic Systems, edited by H. Fukuyama and T. Ando 
(Springer- Verlag, Berlin, 1992); Mesoscopic Phenomena in Solids, edited by B.L. Alt- 
shuler, P.A. Lee and R.A. Webb (North-Holland, Amsterdam, 1991). 

[2] G. Griiner, Density Waves in Solids, (Addison- Wesley, Reading MA, 1994); Charge 
Density Waves in Solids, edited by L. Gor'kov and G. Griiner, (Elseviers Science, Am- 
sterdam, 1989). 

[3] M.I. Visscher and G.E.W. Bauer, Phys. Rev. B 54, 8487 (1996). 
[4] B. Rejaei and G.E.W. Bauer, Phys. Rev. B 54, 8487 (1996). 

[5] H.S.J, van der Zant, O.C. Mantel, C. Dekker, J.E. Mooij, and C. Traeholt, Appl. Phys. 
Lett. 68, 3823 (1996). 

[6] Yu. I Latyshev, B. Pannetier, and P. Monceau, Europ. Phys. Journ. B 3, 421 (1998). 

[7] N. Markovic, M.A.H. Dohmen, and H.S.J, van der Zant, Phys. Rev. Lett. 84, 534 (2000); 
O.C. Mantel, F. Chalin, C. Dekker, H.S.J, van der Zant, Yu.I. Latyshev, P. Pannetier, 
and P. Monceau, Phys. Rev. Lett. 84, 538 (2000). 

[8] Yu. I. Latyshev, O. Laborde, P. Monceau and S. Klaumiinzer, Phys. Rev. Lett. 78, 919 
(1997). 

[9] A. A. Sinchenko, Yu.I. Latyshev, S.G. Zybtsev, and I.G. Gorlova, J. Exp. Theor. Phys. 
86, 1001 (1998) [Zh. Eksp. Teor. Fiz. 113, 1830 (1998)]; A. A. Sinchenko, Yu.I. Latyshev, 
S.G. Zybtsev, I.G. Gorlova, and P. Monceau, Phys. Rev. B 60, 4624 (1999). 

[10] C.W.J. Beenakker, Rev. Mod. Phys. 69, 731 (1997); A.D. Stone, P.A. Mello, K.A. Mut- 
talib, and J.-L. Pichard, in Mesoscopic Phenomena in Solids, edited by B.L. Altshuler, 
P.A. Lee and R.A. Webb (North-Holland, Amsterdam, 1991) 

[11] P.A. Mello and A.D. Stone, Phys. Rev. B 44, 3559 (1991). 

10 



[12] H.S.J, van der Zant, private communication. 

[13] S.N. Artemenko, Sov. Phys. JETP 84, 823 (1997) [Zh. Eksp. Teor. Fiz. Ill, 1494 
(1997)]. 

[14] For simplicity we have taken M to be energy independent. M is then evaluated for 
particles at the Fermi surface. This is a good approximation when the dwelling time of 
the particles on the island is small compared to h/A. 

[15] P.A. Mello, J. Phys. A 23, 4061 (1990). 

[16] O.N. Dorokhov, Phys. Rev. B 37, 10526 (1988). 

[17] PW. Brouwer, Ph.D thesis, Leiden (1997). 

[18] M.I. Visscher and B. Rejaei, Europhys. Lett. 43, 617 (1998). 



11 



FIGURES 

FIG. 1. Schematic pictures of a disordered metal in contact with CDW electrodes in the (a) 
isotropic and (b) anisotropic limits. % denotes the position of the CD wave fronts. 
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